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ABSTRACT 

We revisit the high-energy spectral cutoff originating from the electron-positron pair creation in 
the prompt phase of gamma-ray bursts (GRBs) with numerical and analytical calculations. We show 
that the conventional exponential cutoff should be drastically modified to a steepened power-law in 
practical observations that integrate emissions from different internal shocks. Since the steepening is 
tiny for observations, this "smearing" effect can generally reduce the previous estimates of the Lorentz 
factor of the GRB outflows. We apply our formulation to GRB 080916C, recently detected by the 
LAT detector on the Fermi satellite, and find that the minimum Lorentz factor can be ^ 600 (or 
even smaller values), which is below but consistent with the previous result of ~ 900. Observing the 
steepening energy (so-called "pair-break energy") is crucial to diagnose the Lorentz factor and/or the 
emission site in the future observations. 

Subject headings: gamma rays: bursts — gamma rays: theory — plasmas 



1. INTRODUCTION 

Gamma-ray burst (GRB) is one of the most mysteri- 
ous objects in the universe. The typical energy of its 
prompt emission is ~ 100 — 1000 keV, which corresponds 
to the breaking point of the br oken power-law d escrib- 
ing the energy spectrum well ([Band et al.l Il993t l. Al- 
though there has been many works on the mechanism of 
prompt emission since its discovery, it is still unknown. 
Combined with the optically-thin synchrotron emission 
mechanism, the internal shock model has been most 
wide ly discussed to explain the prompt emission ( see, 
e.g. jMfezaroil2006t[Zhan3l2007tlFMi fc Piranll2008l for 
reviews). In fact, this model can reproduce the behavior 
of complicated light curves (e.g.. iKobava shi et al.! Il997l : 
iDaigne fc Mochkovitchlll998t iBosniak et al.., 2008.1 Col- 
lisions among inhomogeneous outflows lead to the shock 
formation that converts the kinetic energy of the outflow 
to the internal energy. We can interpret the nonthermal 
gamma-ray emission as the emission of electrons which 
are accelerated at the shocks, where inhomogeneous out- 
flows characterized by multiple sub-shells can easily pro- 
duce highly variable light curves. 

High-energy emission of GRBs was already detected 
by EGRET before Fermi, but it has not been clear 
where is the end point of the hig h-ener gy tail (e.g., 
Schneid et all Il992t ISommer et al.ril994 iHurlev et all 



19941 : ISchneid et al.l 1199 5^ If the spectrum extends to 



higher-energy range, sufficiently high-energy photons 
cannot avoid the electron-positron pair creation process. 
Thus it has been expected that there should be an 
cutoff due to this process, although there was no 
observational evidence on it. The pair creation cutoff 
energy generally depends on the emission radius r and 
the bulk Lorentz factor of the outflow F. Hence, we 
can obtain information on these quantities once it is 
observed. Therefore, many authors investigated the 
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possibility to extract the information from the cutoff 
energy fe.g.. iBaring fc Hardind 119971: iLithwick fc Saril 



200H iRazzaaue et all 120041: Murase fc lokal l2008t 



Gupta fc Zhang||2008HGranot et al.ll2008f l. 

However, these studies focused on the high-energy 
emission from one emission zone, neglecting the time- 
evolution of physical quantities during the emission. In 
realistic observations, it is likely that the observed emis- 
sion comes from many emission regions, and physical 
quantities should be time-dependent even during one 
pulse produced by one emission zone. Hence, the cut- 
off behavior in spectra should be affected by both of 
(1) the time-evolution effect during one pulse and (2) 
the time-integration effect over many pulses. On the 
former point (1), the time- variation of properties on 
the photon field during the sub-shell-crossing time scale 
{St/{1 + z) ~ l/c) is significant and comp arable to that 
durin g the total dur ation of the emission (jCranot et al.l 
120081 : iBosnjak et al.l[2008 >l. which can alter the resulting 
high-energy spectrum. The cutoff also evolves reflecting 
the time-evolution of the emissio n site and the Lorentz 
factor of a sub-shell ()Baring||2006D . The latter point (2) is 
also important, since the observed high-energy spectrum 
should be different from that expected in the one-zone 
treatment unless we can separate contribution of each 
emission zone. More specifically, in the internal shock 
model, we would see the superposition of the emission 
from many collisions among different sub-shells during 
the total duration. Especially, if a flux of a single pulse 
is small, we cannot avoid integrating contribution of all 
the pulses to observe the burst. In such a case, it is im- 
portant to study the effect of a lot of internal collisions. 

In this paper, we focus on the above time-integration 
effect over many pulses in the internal shock model, and 
examine whether we can extract information about the 
bulk Lorentz factor of GRBs from high-energy GcV-TeV 
spectra. We model the high-energy emission from the 
multiple emission regions caused by multiple collisions, 
and calculate energy spectra numerically. We need to 
consider a lot of sub-shells with various Lorentz factors 
to explain complicated and irregular light curves. Our 
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numerical approach can treat this situation weU, that 
is vahd even when the variance of the distribution of 
sub-shells' Lorentz factors is large. This is an advantage 
over the analytical approach which is v alid only when 
the variance is small enough (see, e.g., 'Li"2008). In fact, 
the variance should be large in the internal shock model, 
in order that the sub-shells' kinetic energy is efficiently 
converted into the radiation energy a s prompt emission 
(|Kobavashi et al.lll997l : lKobavashi fc Sari 2001). For cal- 
culations, we simplify the emission from each collision 
using a rather phenomenological procedure. This may 
be too simple but enough to study the smearing ef- 
fect by many internal collisions. Other observational ef- 
fects coming from radiation mechanisms will be discussed 
later. 

In ^ we explain the calculational method used in 
this work. We consider multiple sub-shells and integrate 
the energy spectrum of emission which comes from each 
merged shell. In fJH we show the results on the energy 
spectrum of GRBs. In 21 we discuss the possibility 
to extract information about the Lorentz factor of the 
GRB outflow from observations. Moreover, we discuss 
the ca se of GRB080916C, which is an interesting event to 
study (|Zou et al.ll200aiWang et al.ll2009l:[Zhang fc Pe'erl 
\20m . We also compare our numerical results with the 
analytical results. 

2. METHOD 

The inhomogeneous outflow is characterized by mul- 
tiple sub-shells. For simplicity, we only consider com- 
pletely inelastic collisions between shells, and neglect the 
interaction via pressure waves (which is a good approxi- 
mation if the out flow energy is dominated by t he kinetic 
part as shown in iDaigne fc MochkovitchI (|2000D . The ki- 
netic energy of the shells is converted to the radiation en- 
ergy via electron acceleration in the internal shock model, 
where we may phenomenologically give the energy spec- 
trum of emission from each merged shell (see ^2.ip . In 
Appendix, we summarize the dynamics of multiple colli- 
sions w hich was similar to that given in lKobavashi et al.l 
()1997f ). As described in H2.2[ we can evaluate the cut- 
off energy of the spectrum, which originates from the 
electron-positron pair creation process. Hereafter, we ex- 
press physical quantities that are defined in the comoving 
frame of the outflow with primes (e.g., photon energy e') 
and those which are defined in the source (or laboratory) 
frame without primes (e.g., e). 

2.1. Radiation from A Merged Shell 

The internal energy E-mt is determined when two shells 
merge (see Appendix). This energy is divided into the 
energy of protons, electrons and magnetic fields. The pa- 
rameter ep indicates the fraction of the energy in (ther- 
mal and nonthermal) protons. Similarly, the energies of 
(nonthermal) electrons and magnetic fields are charac- 
terized by 60 and ee, respectively. For simplicity, we just 
assume ec — £b — = 1/3. In addition, we consider the 
fast cooling case, where almost all the energy of electrons 
accel erated in each merged s hell is released as radiation 
(e.g. jGuetta fc Granot|[2003l ). 

Observed energy spectra of prompt emiss ion are usu- 
ally fitted well by broken power-law ones (Band et al.l 
[1993,) ■ This broken power-law spectrum is characterized 
by low-energy power-law index a, high-energy power-law 



index (3 and the break energy £b- Therefore, we assume 
that the energy spectrum for each collision is expressed 
by the broken power-law, as expected for the synchrotron 
and/or inverse Compton emission from nonthermal elec- 
trons. Although there may be an additional component 
in the high-energy spectrum (e.g.. [Gonzalez et al...2003l . 
see also later discussions in ^^i^, we approximate the 
spectrum using the broken power law spectrum for the 
demonstrative purpose. We here adopt a = 1, (3 = 2.2 
and £b = 300 keV, which are typical values obtained from 
observations. Given the spectral shape, we can normalize 
a photon distribution function dn/de by the total radia- 
tion energy ec-Eint with 



ec^int = (1 + z)V 



, dn 
dee — , 
de 



(1) 



where V is the volume of each emission region (see Ap- 
pendix). The distribution function of photons in the 
source frame is given by 



dn 
de 




(2) 



Here, K is the normalization factor which is determined 
by Eq. (1). For our calculations, we use the minimum 
energy of 300 eV (which is likely to come from the self- 
absorption process) and the maximum energy of 10^^ eV 
(which roughly corresponds to the most optimistic en- 
ergy of accelerated electrons) before considering photon 
attenuation by pair creation. Note that the results in our 
study are insensitive to these values as long as a < 2 and 
P>2. 

2.2. The Pair-Creation Cutoff 

High energy photons cannot escape from a merged shell 
because of the pair creation processes such as 77 e'^e~ 
and 67 — » ee"'"e~. The most important proce ss is 77 
e+e" under typical conditions of GRBs (Razzaaue et al.l 
l2004f ). We calculate the pair creation cutoff energy for 
high-energy photons by considering only 77 e'^e~ . 

The optical depth r^-y for this process at some energy 
e' can be calculated for a given photon spectrum. We 
assume the power-law spectrum in ij2.1l and calculate the 
optical depth in the comoving f r ame of the outflow as 
follows ("Gould fc Sch^ideJflW . 'Light man fc Zdziarsk: 
1987; Svcnsson 198Jt lLithwic k & Sari 2001: Baring 200^ 
Murase fc Iokall200. 



'77 



(e') mn'{e')aTA', 



e = 



(TOoC^)^ 



(3) 



where the number density of photons whose energies are 
larger than e' is given by 



n'(e') ~ K 



^ ] de'. 



(4) 



Here e' is the energy of a photon which interacts with the 
photon of energy e' at the pair-creation threshold, and crx 
is the Thomson cross-section. Note that equation (U) is 
valid when e' satisfies e' > e^, and we may approximate 
the optical depth as Tj^{e') ~ T-yji^h) for < if 
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spectrum is typical a ~ 1 (see for more discussions in ^ 
and 2]) . The number density of photons in the comoving 
frame n' relates with the number density in the source 
frame n as n'{e') ~ n(£'r)/r because the isotropically 
distributed photons, whose energy is e' in the comoving 
frame, is blueshifted to e ~ Fe' on average. is the 

numerical factor that depends on the power-law index 
and this factor decreases with /3; the values are = 
11/90 ~ 0.12 and = 7/75 - 0.093 for /3 = 2 and 
/3 = 3, respectively ()Svenssonlll987[ ). We can use ~ 



7(/3 - 1)/ [6/35/3 (/3 + 1)] for 1 < < 7 for isotropically 



distributed photons (|Baring||2006[ ). We approximate the 
width of the emission region A' by using the width of 
the merged shell, that is. A' ~ l'^ (see Appendix). The 
pair creation cutoff energy e^^^ is defined as the energy 
at which the optical depth becomes unity, 7'77(£cut) = 1- 

For sufficiently small emission radii or Lorentz factors, 
the Thomson optical depth tt may exceed unity. Even 
in such cases, we always expect that the pair creation 
cutoff should be larger than ~ FmeC^/(l 4- z). How- 
ever, the Compton scattering is ver y important when 
Tt ^ 1 (see, e.g.. iThompsonl Il994t iPe'er fc WaxmanI 
12004 iRees fc M eszarodI I2005D . The photon spectrum 
may have a thermal or quasi-thermal c omponent (see, 
e.g.. IPe'er fc WaxmanI [200I iPe'er et"al.l i2006). It may 
deviate from the powe r-law shape and even have the 
pair annihilation bump (jloka et al.ll2007l : iMurase fc lokal 
|2008() . In this work, we treat the cutoff as Ecut = ^ann = 
VmeC^/{l + z) for T^^(FmeC^/(l + z)) > 1 for sim- 
plicity. Although this may not be a good approxima- 
tion, we expect that this treatment is enough for our 
purpose since we do not have so many collisions with 
r-y^(FmeC^/(l -l- z)) > 1 for our adopted parameter sets. 

For sufficiently large emission radii or Lorentz factors, 
both the Thomson optical depth and the pair creation 
optical depth may be very small, where no pair cre- 
ation cutoff exists because of Tj^{e') < 1 for any en- 
ergy e'. Then, a significant fraction of > TeV pho- 
tons escapes from the source without internal attenu- 
ation. These photons should be attenuated via pair 
creation by the cosmic microwave background (CMB) 
and cosmic infrared background (CIB) photons (see, e.g., 
iMurase et al.ll2007l and references therein ). We include 
this ef fect by using the "low-IR" model of iKneis ke et al.l 
()2004|) . Note that recent observations of TeV blaz ars im- 
ply that th e low-IR model is favored (Aharonian et al.l 
I20M 120071 ). The created pairs can affect the primary 



the CMB photons ( 


Dai fc Lull2002l:lRazzaQue et al.ll2004l: 


IMurase et al.l 12007 


). However, this effect leads to rel- 



atively minor effects during the prompt emission and 
can be important in the late time only when the in- 
tergalactic ma gnetic field is weak enough ([Murase et"an 
[2OO8; Takahas hrerani2008t iMurase et al.ll2008f l. Hence, 
we neglect this possible sub-dominant contribution from 
pair echo emission in this work. 

2.3. Parameters for Shells 

We have the following parameters in our calculations: 
the initial length between two shells d, the initial width 
of a shell I, the number of shells N, the distance from 
GRB to the Earth D (or the rcdshift z), the fraction of 
energy that goes to nonthermal electrons ec, the initial 
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Fig. 1. — The £_Fe spectrum of prompt omission from each 
two-shell collision where seven spectra are randomly selected from 
dozens. The number of shells is set to N = 100 and the distri- 
bution of initial Lorentz factors is log-normal with Fq = 500 and 
an initial fluctuation A = 1, which are defined by equation I)B1|| . 
Other parameters are the initial distance between shells d = 5 X 10^ 
cm, the initial thickness of a shell I = 10^ cm, the initial radius of 
the innermost shell rinj = 10^^ cm, the energy fraction of electrons 
ec = 1/3, and the redshift z = I. We adjust the initial number 
density rip so as to have the mean luminosity of Lj = 10^^ erg 
. The sensitivity curves of LAT on Fermi and MAGIC are also 
shown. 



10=Jerg s"' (with GIB) ■ 
10°;erg s] (wltti GIB) 
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MAGIC 



energy, log(i; [eV]) 

Fig. 2. — The time-integrated spectra for different luminosities. 
We calculate three cases with almost the same cutoff energy but 
different luminosities = 10^^, 10^^ and 10^^ erg s~^. The 
fluctuation of the initial Lorentz factor A is 0.5 for = 10^'^ erg 
s~^, 0.7 for 10^^ erg s~^ and 1.6 for 10^^ erg s^^, respectively. 
Other parameters are the same as figure [T] except for the number 
density np. The number density is determined to set the luminosity 
to the above values. The slope of e"'''^ is shown with the dot-dot 
line for comparison. The sensitivity curves of LAT on Fermi and 
MAGIC are also shown. 

number density of protons in the innermost shells rip, 
and the most inner radius of subshells rjni (see also Ap- 
pendix). We adjust the number density rip so as to have 
the given isotropic luminosity L-y. In addition, we as- 
sume the log-normal distribution of shell's bulk Lorentz 
factors F, whose average value is Fq, with the fluctua- 
tion A. Note that Fq is different from the mathemat- 
ically strict mean value (see Appendix for definition of 
the log- normal distribution). In this paper, the ACDM 
cosmology with fim — 0.3, f^A = 0.7 and Hq — 70 km 
s^^ Mpc^^ is also adopted. 

3. RESULTS 

In Figure (TJ we show the energy spectra of emission 
from each two-shell collision using the method in f}2] We 
only show seven spectra randomly selected from a GRB 
to make the figure easy to see, although there are dozens 
of spectra for one realization. We can see the electron- 
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Fig. 3. — The time-integrated spectra for the different fluctuation of Lorentz factor A = 1, 0.6, 0.2. The luminosity is set to L^/ 
erg in these flgures. Other parameters are as foUows. (a) z = 1, I = 10* cm, = lO'' (b) z = 1, I = 10^ 

500, d 



I = 10* cm, N = lO'^ (d) 



0.3, I = 10^ cm, AT = 10-^. We have To 



5x1, 



10^^ cm and ec 



Iq52 

N = 10^ (c) 'z = 0.3, 
1/3 for aU cases. The 



dot-dot line shows the slope of e " '^ for comparison. We consider the effect of CIB in TeV energy range. The sensitivity curves of LAT on 
Fermi and MAGIC are also shown. 

the effective area under the criterion that at least 5 
photons are collected. The sensitivity curve of MAGIC 
is roughly estimated from the effective area for a zenith 
angle of 20° with the criterion that at least 10 photons 
are collected, al though the actual detectability requires 
careful analysis (| Albert et al.ll2007l ). 

In Figure O the three time-integrated spectra are 
shown whose isotropic luminosities are 10^'^ erg s~^, 10^^ 
erg s~^ and 10^^ erg s~^. The fluctuation of the initial 
Lorentz factor A is 0.5 for the case of 10^"^ erg s~^, 0.7 
for 10^^ erg s^^, and 1.6 for 10^^ erg s~^. We choose the 



index defined between 10 eVand1Q„eV 
ip eVand10^"eV 

0.3 0.4 0.5 0.6 0.7 0.8 0.9 1 
fluctuation of ttie initial Lorentz factor, A 



Fig. 4. — The change of the photon index vs. the fluctua- 
tion of the initial Lorentz factor A. The change of the photon 
index indicates the difference from the initial high-energy power- 
law index f3. The photon index is deflned in two different energy 
ranges, 10* eV - lO^eV and lO^eV - The spectral index 

becomes larger as A becomes larger. The solid and dashed lines 
are quadratic curves obtained from the least square flttings. 

positron pair creation cutoff in each spectrum. One may 
expect to obtain many informations about the GRB if 
we can observe the cutoff energy. However, we have to 
integrate these spectra unless we can observe these spec- 
tra separately due to the lack of photon numbers. In 
figure [1] all the spectra are below the sensitivity curve of 
Fermi (the cutoff energy may be observed by MAGIC. 
But note that MAGIC takes a few dozens of seconds to 
turn the detector to the direction of a GRB after an alert 
from a space telescope.). In this case, only the integrated 
spectrum is observable and the cutoff is smeared out as 
we explain below. 
Here the sensitivity curve of Fermi is calculated from 



values of A so as to set the spectrum steepenings are at 
the same energy (see below). Other parameters are de- 
scribed in the caption of the figure. We plot the spectra 
with and without the CIB effect. 
From Figure [21 we can find following features: 

(1) Most importantly, there is no sharp cutoff that orig- 
inates from the pair production in Figure [21 Instead, the 
spectral slope becomes steep above a few hundred MeV. 
We define the begining point of the steepening as the 
pair-break energy. This is because the steepening arises 
from the superposition of the pair creation cutoff in each 
collision that is plotted in figure [H The exponential cut- 
offs are smeared by the time-integration and result in a 
steepening of the power-law spectrum. Note that this 
feature is observed with other distributions of the initial 
Lorentz factor. 

(2) The steepening effect is conspicuous for the small 
fluctuation of the initial Lorentz factor A (corresponding 
to the high luminosity case = 10^^ erg s~^) while it 

■* Shown at |http: / /www-glast.slac. Stanford. edu/software/IS/glast-| 
lat _ performance.htm. 
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is not for large A. This is because the released energy 
decreases more rapidly for smaller A as the collision ra- 
dius of two-shells becomes larger and accordingly as the 
cutoff energy of a two-shell collision becomes higher. We 
discuss this effect analytically in 21 The steepening ef- 
fect may be observed by Fermi when the luminosity is 
larger than ~ 10^^ erg s~^. 

(3) Although the steepening is small for large A at the 
pair-break energy, the slope gradually becomes softer at 
larger energy (see also figure [3]). This is because the fluc- 
tuation of Lorentz factor evolves smaller with time as 
shells collide with each other, asymptotically approach- 
ing the small A case (see 21 for more discussions). 

(4) The slope of the spectrum becomes hard, return- 
ing back to the original high-energy power-law index (3 
around TeV energy. This means that the spectrum is 
produced by a two-shell collision without the pair cre- 
ation cutoff, i.e., by low-density shells. We can roughly 
estimate the hardening energy (or the transparent en- 
ergy) as ~ r2(mec2)^/(l + z)eh (see Eq. dS])) in the fol- 
lowing way. We assume the photon density is decreasing 
monotonously above the break energy eb- Then, as the 
photon density becomes low with large collision radius 
and/or small internal energy, the low energy range of 
the spectrum below eb becomes relevant for the pair cre- 
ation as target photons. However the photon number is 
insensitive to the photon energy below ei, for the typical 
low-energy power- law index a = 1, and hence the optical 
depth to the pair creation becomes almost constant for 
the photon energy above ~ r^(meC^)^/(l -I- z)eh- This 
is why the slope becomes hard above this energy. Unfor- 
tunately, the hardening effect would not be observed by 
the effect of CIB (see below). 

(5) There is a cutoff around 10^^ eV. This is not the 
pair creation cutoff as discussed above but the cutoff due 
to the maximum energy of electrons that we put for con- 
venience. Note that the maximum energy of electrons 
itself depends on the details of acceleration mechanisms 
that are not in focus of this work. 

(6) If we include the effect of pair-creation interactions 
with CIB, we have a cutoff around TeV energy. This ef- 
fect masks the hardening effect and the maximum energy 
cutoff discussed above. 

In Figure [3l we show the time-integrated spectra for 
different fluctuation of initial Lorentz factor A from 1 to 



0.2 with the same luminosity Lj — 10 erg s 



We 



consider the effect of CIB in TeV energy range. The 
spectrum becomes steep above a pair-break energy. A 
pair-break energy is low for large A because the collision 
radius is small and hence the cutoff energy is small for 
each collision. The steepening is observable for large A 
even if a GRB occurs at high z ^ 1 ((a) and (b)). We 
can observe the steepening for small A if the GRB occurs 
at low z ((c) and (d)). We also compare the cases of 
N = 103 {I = lO^cm) with that of iV = 10^ {I = lO^cm). 
The shells collide with each other in the inner radius for 
= 10^ than for A^ = 10^ because the separation is 
small (The separation between shells are small when the 
thickness of the shell is small. This is because we assume 
the separation and the thickness are comparable.). It 
leads to the smaller pair creation cutoff in each collision 
and the pair- break energy becomes smaller. 

In figure IH we show a change of the power-law index 
from (3 — 2.2 due to the effect of the smearing for the 




0.1 

pulse interval, At [sec] 

Fig. 5. — Pulse interval distribution. The parameters of the 
initial Lorentz factor distribution are Fq = 500 and A = 1. Other 
parameters are the same as the calculation of figure [Sfb). We fit 
the data of pulse intervals by the log-normal distribution with the 
least square method. 




9 10 11 

pair-break energy, log(Epb[eV]) 

Fig. 6. — The pair-break energy vs. standard deviation of the 
pulse interval distribution. The pair-break energy is defined by 
the energy where the spectrum becomes steeper than /3 by 0.05. 
We calculate these values by changing Fq and A, which are the 
parameters of the initial Lorentz factor distribution. The solid line 
is written for Fq = 300 changing A from 1 to 0.2. The dashed line 
is written for Fq = 500 changing A from 1 to 0.2. The dotted line is 
written for Fo = 700 changing A from 1.3 to 0.4. The dot-dashed 
line is written for Fq = 900 changing A from 1.5 to 0.6. Other 
parameters are the same as that in figure [Sfb) 

case of figure [2Ja) (the slope of figure [3]^c) is same with 
that of figure [3l^a)). We calculate the smeared power-law 
index by two ways. One is defined by the power-law in- 
dex between 100 MeV and 1 GeV with the least square 
fitting. This energy range corresponds to the observable 
energy band in figure ^[&) . Another is defined by the 
energy range between 1 GeV and 10 GeV, which cor- 
responds to the observable energy range in figure [3I^c). 
Figure |4] shows that the slope in a fixed energy region is 
steeper for larger A. The reason is that the pair-break 
energy is lower for larger A (see the previous paragraph) 
and the slope becomes softer as the fixed energy range 
is more separated from the pair-break energy (see the 
feature (3) obtained from figured]), although the steep- 
ening is smaller at the pair-break energy for larger A (see 
the feature (2) obtained from figure [2]). We have more 
discussions in 21 

We also calculate pulse intervals. It is known that 
the distribution of pulse intervals can be fitted by a log- 
normal distribut ion although there is a small excess at 
long in tervals (McBrecn et al. 1994: Li fc Feni morelll996l : 
iNakar fc Piran 2002; loka fc Nakamura ,2002 ). We cal- 
culate pulse intervals for various average Lorentz fac- 
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Epi,: pair-break energy 



Fig. 7. — Schematic picture of the GRB spectrum. The cutoff 
is smeared by the time-integration although each energy spectrum 
for two-shell collision has a cutoff. There is no cutoff in the time- 
integrated spectrum but the slope becomes steep above the pair- 
break energy. This steep slope (the smeared cutoff) is made by 
the summation of different luminosity emission. The pair-break 
energy is almost same as the cutoff energy which originates from a 
two shell collision at the most inner radius (ecuti in this picture). 

tors Tq and fluctuation amplitudes A (see Appendix for 
the definitions). We plot the result for Fq = 500 and 
^ = 1 in figure [5] as the histogram of the logarith- 
mic time intervals. The result is averaged for 100 dif- 
ferent random realizations. A pulse interval is defined 
by an interval between each begining time of a pulse. 
We fit the histogram by the log-normal distribution us- 
ing the method of least squares. A log-normal distri- 
bution is characterized by ^ and as f{At)d\n At = 
exp(-(ln(At) - ^)^/2a^)/V2^d\nAt. 

In figure [SI we plot the correlation between the pair- 
break energy £pb and a. We define the pair-break energy 
by the energy where the power-law index becomes steeper 
than P by 0.05. We calculate Epb and a for various Tq 
and A fixing the luminosity. In this figure, Cpb becomes 
smaller for larger A. The cutoff energy of two-shell colli- 
sion is smaller when the fluctuation is large because the 
collision radius is smaller and the density of photons is 
larger. It means £pb becomes small when A is large. It is 
difficult to observe £pb when the fluctuation is small by 
the effect of CIB. We discuss the possibility to extract 
the information from this plot in 21 

4. DfSCUSSlON 
4.1. On Estimate of Lorentz factors 

Even though a pair creation cutoff energy may give us 
useful information as we mentioned in ^ it does not 
seem easy to obtain information when the cutoff is hid- 
den by the smearing effect (see figure [7]). If our model 
assumptions are valid, it is difficult to see the obvious 
cutoff feature in the internal shock model. A clear cutoff 
feature should be in the very high energies of ^ GeV, 
but it is likely to be hidden due to attenuation by the 
CIB. Then, can we probe the Lorentz factor from the 
GRB spectrum? We think that we can still extract use- 
ful information on it. 

One possibility is that we may constrain the Lorentz 
factor from the absence of the pair creation cutoff. Since 



each cutoff is produced by the target photons in each 
two-shell collision, we may estimate each spectrum by 
dividing the time-integrated spectrum by the pulse num- 
ber of the light curve. However, this method usually 
overestimate the target photon density because several 
pulses sometimes overlap. In addition, observed pulse 
durations and intervals may not be ideal ones (that are 
directly calculated in the internal shock model) due to 
the overlap. In particular, high-energy gamma rays re- 
lated to the maximum cutoff come from the outer colli- 
sion with a wide pulse, which would be buried in many 
other spiky pulses. Hence, we expect that the cutoff 
smearing effect can generally reduce the conventional es- 
timate of the Lorentz factor. 

One may see the pair-break energy in high-energy spec- 
tra. The pair-break energy is almost the same as the 
minimum cutoff energy in two-shell collisions. Hence, the 
pair-break energy would be useful instead of the original 
cutoff energy, if it is clearly observed. Although its ob- 
servational identification may not be easy, let us consider 
the case that it is accomplished. Then, we can apply 
the c onventional procedure described by various authors 
fe.g.- lLithwick fc Sarill200lt |M urase fc Iokall20M l. From 
equations ((!]) and (0]), we have (for /3 > 2) 



(3-2' 



n'{e'): 



K {\ + z)£b f 

't 13-1 [T^ 



1-/3 



(5) 



(6) 



By substituting these equations into equation ^ , we ob- 
tain the Lorentz factor as 

^ /epbeb(l + 2)2\^ 



V 167rc2eb5i2(^_ 1) j ' ^ > 

where we use V ~ Airr'^A and r 2r'^A T^cSt/{l + z) 
expected in the internal shock model (and we can also ex- 
pect r ~ 2T^cSt/ (1 -f z) if a pulse duration is determined 
by the pulse decay time due to the curvature effect). For 
typical parameters, we have 



r~3.2 X 10^ 



500MeV 



(1 + ^)^ 7777 



St 

OAs 



\ 1053, 



erg 



(/3 = 2.2,eb 300keV), (8) 



where we note that ecE-mt is the released radiation energy 
per collision. 

How could we estimate the Lorentz factor from the 
time-integrated spectrum? First, we could roughly esti- 
mate the released energy from each collision, ec-Eint by 
dividing the time-integrated spectrum by the pulse num- 
ber of the light curve, which is observed in lower energy 
range near Sh, once D is determined by other observa- 
tions. We could also estimate St ^ l/c ^ d/c ^ At ii 
we can know the pulse width in light curves, although it 
is difficult to evaluate these values accurately from ob- 
servations. Then, we could roughly estimate the Lorentz 
factor F from the pair- break energy £pb. One advantage 
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of using the pair-break energy is that the minimum cut- 
off energy, i.e., the pair-break energy, is predominantly 
produced by the inner coUision with a short pulse width, 
which may be directly identified with reasonable values of 
the observed shortest variability time scale. Another ad- 
vantage is the pair-break energy is smaller than the cutoff 
energy in general and it can be smaller than the CIB at- 
tenuation energy. Several pulses still often overlap, which 
may lead to our overestimating the released energy from 
each collision ec-Eint via overestimating the duration and 
interval of pulses from the most inner radii. Hence, as 
long as we can choose a reasonable shortest value of 6t, 
this method could give an upper limit of the Lorentz fac- 
tor (i.e., the true value of the Lorentz factor Ftrue, which 
is obtained considering hidden pulses, is related with the 
Lorentz factor calculated by our formulation Fuppor and a 
conventional estimate of the Lorentz factor Fco as follows, 
Ltruc < Fuppor < Fco). Again, notice that conventional 
estimate generally lead to overestimating the minimum 
Lorentz factor when we use short time scales of St since 
true pulses related to the high-energy emission may be 
wide and hidden by other pulses. 

In the internal shock model, the Lorentz factor is likely 
to have a large dispersion, i.e., large A. Since the disper- 
sion of the Lorentz factor A is reflected by the dispersion 
of the pulse interval cr, we may expect some correlation 
between a and Sph (see figure [6]) . Although this a is 
not the observed dispersion itself which depends on al- 
gorith ms determining pulse widths (e.g.. iLi fc Fenimord 
it suggests that we can obtain some potential 
clue to the average value of the Lorentz factor from high- 
energy gamma-ray observations. 

4.2. Comparison with Analytical Calculation 

We can calculate a power-law index produced by 
the smeared cutoff analytically if an initial variance of 
the sub-shells' Lorentz factor distribution (Tp q satisfies 

Up Q < Fq. Here Fq is an average Lorentz factor and 
(Tp g is defined in the source frame. Other assump- 
tions are the same as the numerical calculation in this 
study. In this case, the variance of the shells' veloc- 
ity distribution tr^, which is defined in the center of 
mass, evolves with a radius of the expanding shell r as 
ay oc r~^^^. The internal energy that orig inates from 
two-shell collision is proportional to fjBeloborodovl 
I2000t) and this energy decreases as E-mt oc oc r ^'^ 
with radius (jLi fc Waxm an 20081. The internal energy 
is converted to nonthermal energy of electrons that radi- 
ate photons. We discuss a situation that a photon whose 
energy is Ecut annihilates with a photon whose energy is 
e (> £(,). Here, e is the typical energy of target pho- 
tons. Then the target photon number density is writ- 
ten as n oc Bints'^- '^/ATTr^AT'^ oc r'^/^i'^-'^ A''^ . We 
can calculate how the optical depth depends on r and e, 
T^j{e) oc narA oc £i-/37--8/3. The target photon energy 
scales as e oc r^/'^^^^f^). High energy photons with en- 
ergy e that satisfies e > (rnieC^)'^ / {1 + z)e produce pairs 
in interaction with photons. We evaluate Ecut using the 
minimum energy for the pair-creation threshold. Then, 
the cutoff energy scales as Ecut (x r~*/'^^^~'^) The inter- 
nal energy which associates with the fluctuation of shell 
velocities decreases with shell's radius as Ei^t oc r~^/^ 



()Li fc WaxmanI [20081 ). Hence Eint relates with the cut- 
off energy as Eint oc e^]^^''^'^. The energy flux at Scut is 
in proportion to E^t and e^ut if assume the energy 
spectrum that originates from two-shell collision satisfies 
the power-law and the break energy (or peak energy) of 
each shell's spectrum does not change. The high energy 
part of time-integrated spectrum would be 

£i^e(£) OC -Binte^^'^ oc £— . (9) 

The pair- break photon index /3pb (see figure [7]) becomes 
(5/?-l)/4 and the change of the index is (/3-l)/4 (= 0.3 
for f3 = 2.2). 

For large fluctuation of the Lorentz factor Cp q > Fq, 
i.e., A > I, we cannot apply the above analytical results. 
Actually we have shown that the steepening is more mild 
for larger A in fJSl 

In figure [3l we plot the time-integrated spectrum. The 
begining energy of steepening depends on the fluctuation 
of the initial Lorentz factor A. The steepening begins at 
low energy when A is large. The collision radius is small 
for large A and the density of photons becomes large. 
This leads to the small cutoff energy of two-shell colli- 
sion and the steepening begins at lower energy. For large 
A, the steepening is small at the pair-break energy (see 
fJ3]), but the slope becomes steep gradually at high en- 
ergy because the fluctuation of Lorentz factor becomes 
small as shells collide with each other and the slope ap- 
proaches the behavior of the analytical solution. In fact, 
i^int follows the analytical solution, oc r"^/^, after some 
collisions. However, there is little difference between the 
numerical calculation and the analytical one even if the 
slope becomes steep. This can be understood as follows 
(see also Photons whose energy is lower than Sh, 

which is deflned in ^2.1[ is important when the cutoff 
energy is large. The number of these photons are small 
to interact with other photons and the cutoff energy be- 
comes large suddenly. In our numerical calculation, we 
approximate there is no cutoff when the cutoff take such 
a large value. On the other hand, we consider a sin- 
gle power-law without the spectrum break at Eh in the 
analytical calculation. In this case, there are many low- 
energy photons and the cutoff energy increases following 
the power-law of the collision radius. We checked the 
numerical calculation gives the same result as the ana- 
lytical calculation in a low A regime when we assume 
the photon spectrum satisfles a single power-law for a 
test numerical calculation. We can understand the be- 
havior of the smeared cutoff in figure [2] similarly. The 
pair- break photon index which is defined between IGeV 
and lOGeV is 0.4 for = lO^^ gj.g g-i^ q.S for L-^ = 10^^ 

erg s~^ and 0.3 for ~ 10^^ erg s^^, respectively. 

4.3. Comparison with Recent Observations of Fermi 

We can estimate the Lorentz factor of a GRB if we 
observe a cutoff in a single pulse spectrum. Then 
one may think that the maximum energy of the ob- 
served photon gives a lower limit of the Lorentz 
factor (iBaring fc HardiMi Il997t iLithwick fc Saril l2001l : 
iRazzaaue et al.l 12004 ). However, as we have discussed 
in MA\ the cutoff smearing effect generally allows lower 
Lorentz factor than the conventional estimates. 

In the case of GRB080916C, 3 GeV photon was de- 
tected by Fermi LAT in first 10 seconds and 13GeV 



8 



Aoi et al. 



photon was detected in 100 seconds (jAbdo et al.ll2009D . 
We can constrain the Lorentz factor using the maximum 
energy of a detected photon. The minimum Lorentz 
factor is 890 when we use fohowing parameters, £b — 
1.17 X lO^eV, e = 3 X lO^eV, f3 = 2.21, z = 4.35, 
Eint = 4.2 X lO^^'erg and St = 2s. Here, we regard the 
St — 2s pulse as a single pulse with energy Ei^t- However, 
there may be a pair-break energy because the spectral 
steepening (by less than (/3 — l)/4 ~ 0.3) is too smah to 
be detected with the current high energy data. When we 
use the high t emporal variabihty ft ~ 100 ms detected by 
INTEGRAL (jCxreiner et al.l[20(M) and E^nt = 4x 10^^ gj-g 
accordingly, the conventional estimate on the minimum 
Lorentz factor gives ~ 1100. However, this value may be 
overestimated. Here, let us apply our estimate using the 
pair-break energy. Then, we can obtain the Lorentz fac- 
tor as - 600(epb/100MeV)^/^'^ for Eint = 7x10^^ gj.g ^nd 
St = 100 ms. The actual Lorentz factor could be even 
lower as long as our choice of St is reasonable, because 
the spectrum is the summation of many pulses' spectra 
in the cutoff smearing picture leading to overestimating 
the single pulse energy Ei^t (see equations ([7]) and (O). 

A hard lag is observed in the observation of 
GRB080916C. It might suggest that the time- 
development of a cutoff energy. Generally, the cutoff 
energy from a inner radius is smaller than that from an 
outer radius. The hard-lag might be explained by mul- 
tiple collision of shells because a photon from a inner 
radius seems to be detected first. However, we did not 
find the time-lag in the different energy range in our cal- 
culation. Photon from a inner radius can be detected 
after a photon from an outer radius because the shell 
is distributed with a moderate width. It is difficult to 
reproduce the time-lag in the current model. 

4.4. Effects of Radiation and Acceleration Mechanisms 

In this study, we do not consider an additional 
high-energy component in the prompt emission spec- 
trum. Such an additional component has been often ex- 
pected in the optically thin synchrotron scenario (e.g., 
iMeszaro^ [2006: Zhang 2007). Theoretically, there are 
the two main classes as high-energy emission mecha- 
nisms, i.e., leptonic and hadronic mechanisms. The 
leptonic mechanisms include synchrotron self-Compton 
(SSC) emission and external inverse-Compton emis- 
sion, which are the most discussed scenarios. High- 
energy SSC emission is produced by relativistic elec- 
trons that radiate seed synchrotron photons theniselves 
(e.g..lPapathanassiou fc MeszaroslllQQl iDai fc Lull2002t 
iGuetta fc Granotll2003ri . In addition, there may be some 
possibilities for external inverse-Compton emission via, 
e.g., upscattering thermal photons in the cocoon. The 
hadronic mechanisms include synchrotron radiation of 
high-energy baryons, synchrotron radiation of the sec- 
ondary leptons generated in photohadronic interactions, 
as well as the photons directly produced from tt*^ de- 
cays. In order to see the baryon synchrotron radia- 
tion, sufficie ntly strong magnetic fields are typically re - 
quired fe.g.. iGupta fc Zhang! [20071 : iMurase et al.ll2008f ). 
Otherwise, photohadronic components would dominate 
over the baryon synchrotron co mponent as long as the 
photon density is hi gh enough (Der mer fc Atovanl[2003 : 
lAsano fc Inouel[2007.) . Hadronic gamma rays can be ob- 
served only when the nonthermal baryon loading is large 



enou gh (e.g., iMurase fc Nagatakil l2006t lAsano fc Inoud 
|2007( ). The above mechanisms can lead to an additional 
high-energy component that can be observed by Fermi, 
and high-energy spectra may differ from those demon- 
strated in this work. Although they are potentially im- 
portant, we postpone detailed discussions in this work 
for our demonstrative purpose. 

Another uncertainty raises from the origin of the cut- 
off in high-energy spectra. For example, the maximum 
energy of generated photons is at most the maximum 
energy of electrons, which is often determined by com- 
paring between the synchrotron cooling time and accel- 
eration time in the optically thin synchrotron scenario. 
If the acceleration efficiency is small enough, the cutoff 
would be determined by the electron maximum energy 
rather than the pair creation process. In addition, other 
processes suc h as the Thomson scatt ering may also be 
relevant fe.g.. iPe'er fc WaxmanI [200^ . Also, the accel- 
eration mechanism affects the electron spectrum itself 
as well as the maximum energy of electrons. Generally 
speaking, the relativistic diffusive shock acceleration sen- 
sitively depends on properties of magnetic fields around 
shocks, so that the spectral index may not be univer- 
sal and th e spectrum may differ f r om the simple power 
law (e.g., l Achterberg et al.' 2001: Blasi fc Victri 20051; 
Lemoine fc Rc venu 2006 ; Nicm icc fc Ostrowski 200^ 
Morlino et all l2007a| bl: lAoi et all 120081 ). These effects 
from radiation and acceleration mechanisms may signifi- 
cantly affect our results of high-energy spectra made by 
the smeared pair creation cutoff. 

5. SUMMARY AND CONCLUSION 

We have investigated the spectral features of the 
electron-positron pair creation in the GRB prompt emis- 
sion to probe the Lorentz factor of the GRB outflow. We 
have studied the time-integrated spectrum from many in- 
ternal shocks with numerical and analytical calculations. 
This is because almost all the authors studied the pair- 
creation cutoff energy expected from one emission region, 
although it could be difficult to observe the emission from 
a single collision. Although the cutoff energy may be ob- 
served for a sim ple pulse from bright bursts and it is 
very useful (e.g., IMurase fc lokal 120081: iGupta fc Zhand 
|2008| ). it is also important to study whether the cutoff 
can be observed or not for bursts that have many pulses 
in their light curves. Fermi would allow us to obtain only 
the time-integrated spectrum in many cases due to the 
difficulty to collect many high-energy photons and the 
overlapping of pulses. 

We demonstrate that the clear cutoff, which comes 
from electron-positron pair creation processes, is hard 
to observe in the internal shock model, when bursts con- 
sist of many sub-shell collisions. Instead, the slope of 
the spectrum becomes steeper above the energy, so-called 
the pair-break energy Sph, which also originates from the 
electron-positron pair creation. We have demonstrated 
how the cutoff is smeared, and compared our numerical 
results with analytical consideration. Possibly, we may 
expect that Fermi observes this pair-break energy in the 
high-energy range. The pair-break energy may be ob- 
servable if the luminosity is greater than ~ 10^^ erg 
and the fluctuation of the initial Lorentz factor A is large 
for z — 1. However, the change of the spectrum slope is 
small in general and it may not be easy to observe the 
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pair-break energy. Nearby GRBs may give us a chance 
to detect the pair-break energy. If it is detected, we may 
have some information on the Lorentz factor of the shell 
emitting photons with ^ epb- 

The cutoff smearing effect can generally reduce the 
conventional estimate of the Lorentz factor. We dis- 
cuss the possibility that the maximum energy of detected 
photons does not provide the best indicator of the min- 
imum Lorentz factor. The pair-break energy is also eas- 
ily missed because the steepening of the spectral index 
is relatively small (less than (/3 — l)/4 ~ 0.3). Then 
the current observations still allow smaller Lorentz factor 
than the previous estimates. We have applied the cut- 
off smearing effect to GRB080916C and suggested that 
the observations are consistent with the Lorentz factor 
of ~ 600 (and even smaller value) that is slightly below 
the previous result of ~ 900. 

The Lorentz factor of the outflow is likely to have a 
large dispersion in the internal shock model. When A is 
large, the resulting spectrum has lower pair-break energy 
where the smearing starts, and vice versa. This suggests 
that dispersion of observed pulses may have some clues 
to estimate the Lorentz factor of the outflow. 

There are differences between the results of the numer- 
ical calculation and those of an analytical one although 
the analytical calculation is a good approximation for 
small fluctuation of an initial Lorentz factor. The steep- 
ening is small in the low energies for large A, compared 
with the analytical calculation. However, numerical val- 
ues approach analytical values in the high energies by 
the effect of multiple collisions. 

Although we have demonstrated that the smeared cut- 
off behavior is characterized by steepened power-law 
spectra, there are a lot of uncertainties in predicted high- 
energy spectra even in the internal shock model. (1) We 
have fixed eb, £c and /3, but changing either of them for 
each collision easily changes the resulting spectra. Even 
the power-law approximation for accelerated electrons 
may not be good, because it depends on the unknown 
acceleration mechanism. (2) There may be an additional 
high-energy component coming from e.g., the inverse- 
Compton scattering and/or baryon- initiated electromag- 
netic cascade processes. (3) The cutoff may originate 
from other processes such as the maximum energy of ac- 
celerated electrons. There is also some uncertainty in 
the high energies where attenuation by the CIB is im- 
portant due to the ambiguity of the CIB. Although the 
smeared cutoff spectrum may easily change if we con- 
sider the above effects, it is commonly expected that the 
obvious cutoff cannot be seen for time-integrated spec- 
tra of bursts with many pulses in light curves, as long 
as we believe the internal shock model. Such a smear- 
ing effect is also expected for other high-energy features. 
For example, a part of the authors provided the recipe 
to diagnose the mechanism of prompt emission for indi- 
vidual pulses and d iscussed the pair annihilation bump 
(jMurase fc loka' 2008). We can expect the smearing ef- 
fect if photospheric emission is also caused by the inter- 



nal shock dissipation. Even high-ene rgy neut r ino sp ectra 
(see for high-energy neutrinos, e.g., iMurasd (|2007| ) and 
references there in) can be modified since the magnetic 
field will be weaker at larger e mission radii that can affec t 
some of the previous results (|Asano fc Nagatakil l2006Y 
The recently launched Fermi can be useful for testing 
the smeared cutoff spectrum. However, it seems that we 
still need many > GeV photons for further investigations. 
Therefore, not only Fermi but also current Cherenkov 
detectors such as MAGIC or planned low-energy thresh- 
old Cherenkov detectors such as 5@5 will be important. 
Once we see high-energy events by these larger area tele- 
scopes, good photon statistics will allow us to study the 
high-energy behavior of prompt emission well. 

Recently, Zhang & Pe'er argued that the emission ra- 
dius of GRB080916C should be larger than r > 10^^ cm. 
Our estimate based on the internal shock model sug- 
gests r > 10^^'^ cm. Since the true Lorentz factor 
can be smaller than our estimate, even smaller emis- 
sion radii are possible. They also argued that the pure 
fireball model may not avoid the photospheric emis- 
sion. However, if the radiation energy has been con- 
verted into the kinetic energy at tt, thermal emission 
may be reduced more. Even if the outflow is initially 
Poynting-dominated, the internal shock model is still 
one of the viable models as long as a significant frac- 
tion of the ma gnetic energy can be convert ed into the 
kinetic energy (|Vlahakis fc K5umlniglll2003D . Recently, 
models that GRB emission is powered by dissipation of 
the Poynting flux energy within the outflow have been 
discussed more and more (jLvutikov fc Blandford 200|; 
iNaravan fc Kumai1l2009f ). Future observations of broad- 
band spectra that are well time-resolved will give us very 
crucial clues to the realistic GRB prompt emission model. 
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APPENDIX 



In this appendix, we explain the dynamics of shells which approximate the inhomogeneous outflows of GRBs. At 
flrst, we explain the two-shell collision. Next, we explain the dynamics of multiple shells. 
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TWO-SHELL INTERACTION 

We consider the situation that a rapid shell collides with a slower shell. We assume two shells merge after collision. 
This assumption is valid when the merged shell becomes cool immediately (i.e. the cooling time by emission is shorter 
than the dynamical time in which a shock wave crosses a shell). We calculate the physical quantities of the merged shell 
defined in the source frame, Lorentz factor r„i, the internal energy E-mt, number density of protons np_m, thickness of 
the shell and area S„i ~ 47rr^ using initial values of Lorentz factor T, density of protons Up, thickness I and area 
E ~ 47rr^ (i.e. volume of the shell is 1^ ~ EZ) of two shells. Physical quantities of a rapid (slower) shell is denoted by 
the subscript r (s). The kinetic energy of two shells converts to the internal energy when two shells collide and shock 
waves occur. Using conservation of momentum and energy, we can calculate the velocity and the internal energy of 
the merged shell. 



m„rs 



(Al) 



Eint = m,c\T, - r„) + msc'(r. - r,„), (A2) 

where m = rUpUpV (rrip: proton mass). This internal energy is radiated immediately. We calculate the energy 
spectrum of the emission in §2.11 

Two shock waves occur when shells merge, a forward shock and a reverse shock. The velocity of a forward (re- 
verse) shock Ffe (Trs) is written as Tf, ~ r™^(l + 2r„/r3)/(2 + r,„/r,) (r,, ~ T^^{1 + 2T,^/r,)/{2 + T^/T,)'^ 

(|Sari fc Piran|[T995t iKobavashi et al.lll997t ). These shocks compress the initial shells and the thickness of the merged 
shell Im is given by: 

l^ = l^^^+l^^^L^. (A3) 

Pis - Ps Py - Prs 

For simplicity, we assume the density of the merged shell becomes homogeneous although there is the contact discon- 
tinuity in fact. We use the averaged density written as: 

rip.m = — -J — ■ (A4) 

EMISSION FROM MULTIPLE SHELLS 

We consider the multiple-shells colliding with each other and the emission from the merged shell. We can calculate 
the energy spectrum of the emission by applying the discussion of the previous section to eac h collision. 

We calculate the dynamics of shells in one dimension numerically (|Kobavashi et al.lll997f ). We consider N shells 
which are labeled by an index i (i = 1, . . . , TV), where the inner shell is labeled by the larger number. These shells are 
characterized by four variables, Lorentz factor Ti, density n^, thickness of a shell h and the initial position of (inner 
part of ) the shell r^. We determine the initial position so as to the length between two shells di{= Vi ~ r^+i — U+i) 
becomes equal (i.e. di is same for all shells). We also assume di is comparable to U and set di — 5 x li. Lorentz 
factors are highly-relativistic value and distributed randomly following the log-normal distribution. The log-normal 
distribution is defined by Tq and the amplitude of the fluctuation A as follows, 

F - 1 e-^'/2 

In = Ax, P(x)dx = — ^^dx. (Bl) 

Lo-l 72^ 



Its mean value is exp(ln(Fo — 1) + and its variance is (exp(A^) — 1) exp(2 ln(ro — 1) + A^){= tXp). At A < 1, we 

have (F — Fo)/Fo ~ Ax, then Fq becomes the mean value and cr/Fo ~ A. At A > 1, we are in the high amplitude 
regime. The specific values of parameters are described in the caption of the figures. Note that we chose the width 
of the shell in order to set the variance time at about 10 or 100 ms. At t = the shells are at the initial positions. 
These shells expand spherically with highly-relativistic speed and density decreases in proportion to r~^, where r is 
the radius of the expanding shell. Strictly, radial velocity spread causes a gradual spread of radial width of shells at 
large radius. Although the shell spreading would be important at large radii, we does not consider the shell spreading 
in this paper for ease of understanding. A rapid shell catches up with a slower one and collides with each other. We 
assume two shells merges after collision and the shell becomes cold immediately via emission. We can calculate the 
internal energy which occurs in the jth collision i?int,j and the energy spectrum of photons dnj/de (see ^A^- We also 
calculate the position Rj and the time tj. An observer at D away from the central source will begin to detect the 
emission at a time ij^ob — [D — Rj]/ c + tj. 

Shells expand and collide with other shell one after another. This process continues until there is no rapid shell 
behind the slower shell or all shells merge into one shell. Note that we do not consider the interstellar medium (ISM). 
Most shells collide at the small radius where the effect of the ISM is negligible when we consider the large value of A. 
The shell collide at the large radius and the effect of the ISM is important and the shell would be decelerated when 
we consider low A. 



Can we probe the Lorentz factor of gamma-ray bursts? 
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We have the time-integrated energy spectrum drisum/de — Ti-^-^^drij/de^ where A^coii is the number of coUisions. We 
focus on the time-integrated spectrum in this study. We define the duration T of the GRB by the time when we detect 
95 % of the total energy. We define the time-averaged himinosity by L-^ = (1 + z)" [J dee {dnsnra/ de)\ /T. 
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